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ABSTRACT 

We present observations of high-amplitude rapid (2 s) variability toward two bright, compact ex- 
tragalactic radio sources out of several hundred of the brightest radio sources in one of the 30° x 30° 

MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After 
rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for 
this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection 
propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty 
Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to > 10^ km. 

The serendipitous night-time nature of these detections illustrates the new regime that the MWA has 
opened for IPS studies with sensitive night-time, wide-held, low-frequency observations. This regime 
complements traditional dedicated strategies for observing IPS and can be utilized in real-time to fa¬ 
cilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact 
(arcsec) structures in low-frequency radio sources despite the 2' resolution of the array. 

Subject headings: scattering — techniques: interferometric — Sun: coronal mass ejections (GMEs) — 

Sun: heliosphere — radio continuum: galaxies 


1. INTRODUCTION 

Interplanetary scintillation (IPS) — the rapid (<ls- 
10 s) hux density variability of distant compact radio 
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sources due to radio-wave propagation through inhomo- 
geneities in th e ionized solar win d — was discov e red by 


Clarke (1964) and published by I He wish et al. (1964); 
the discovery led to high-time-resolntion instrum entation 


that, in tu rn, led to the discovery of pulsars (Hewish 


et al. 


1968). Interplanetary scintillation is an important 


technique m monitor ing the structure and evolution of 


the solar wind (e.g., Coles 1978| Manoharan & Anan- 


thakrishnan||1990t [jackson et al.||1998p, 

probe of major perturbations caused b; 




|), particularly as a 
l y solar flares and 


2003 

Manoharanl 


), and in constraining the nature 

of background rad 

[io sources (Little X Hewish 

1968). 


frequencies, < 400 MHz, where the effects of plasma in¬ 
homogeneities are most apparent. At fac ilities such as 
the Ooty radio telescope (|Manohara n|2010|), the Mexican 
Array Radio Telescope (MEXART'; Ivlejia-Ambriz et al 
|2010|, and the Solar Wind Imaging hacility, Tbyokawa 
IFT; 


Tokumaru et al. 


2 01 1|) /Solar-Terrestrial Envi - 
ronment Laboratory (ST'EL[ Kojima X Kakinuma||1990), 
observations measure the modulation of hundreds of pre- 
selected radio sources each day and follow solar-wind 
events identified elsewhere. The new generation of low- 
frequency, wide-field radio tel escopes such as the Murchi¬ 
son Widefield Array (MWA; |Lonsdale et al. 2009| |Tin- 
gay et al.|2013|) and Low Frequency Array (LOb AR; 


Haarlem et al 


van 


20^ can undertake both dedicated ob¬ 


servations of IPs lUberoi fc Benkevitch|2()lQ Bisi et al 


2011 Fallows et ah ^13| ) — often at high time resolu¬ 


tion — and make serendipitous discoveries, as discussed 
here. These ground-based radio IPS facilities comple¬ 
ment space-based heliospheric imagers such as STEREO 
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(lEyles et al.|2QQ9l), SD O (jSchou et al.|2012|), and Hinode 
(Culhane et al.’2007), and in-situ solar wind measure- 
ments from other spacecraft. 

The MWA is a low frequency (80-300 MHz) interfer¬ 
ometer located in Western Australia, with solar, helio- 
spheri c, and ionospheric s tudies as one of its major fo¬ 
cuses (Bowman et al. 2013|). The advantages of the MWA 
for IPS are: its southern hemisphere location (all major 
IPS facilities are in the northern hemisphere); its very 
wide field of view (over 600deg^ at 150MHz); its high 
sensitivity for short integrations (^100 mJy/beam rms 
for a Is integration); and its capability to make sub¬ 
second flux density measurements. Here we present the 
first serendipitous, exploratory observations of IPS us¬ 
ing the MWA in its normal imaging mode with 2 s time 
resolution. 

2. OBSERVATIONS AND DATA ANALYSIS 

Data for this investigation were from commensal MWA 
proposalG0009 (“Epoch of Reionisation”) and G0005 
(“Search tor Variable and Transient Sources in the EoR 
Eields with the MWA”). The observations are of one of 
the MWA Epoch of Reionisation (EoR) fields (designated 
EoRO, and centered at J2000 RA=0^, Dec=—27°; see 
Eigure]^ and occurred between 11:30 and 14:14 UT on 
2014 November 06, while the field was within approxi 
mately two hours of transit. The observations (see 


et al. 112015^ for more details) used 40 kHz fine channels 


Loi 


2 s integrations; the correlated data were written to files 
of 2 min each, with gaps of 16-24 s between adja cent files. 

Images were searched for East Ra dio Bursts (jLorimer 
let al.|]2QQ^ |Thornton et al.||2013|) using a processing 
pipeline which will be described elsewhere (Tingay et ah, 
in preparation). A by-product of this processing is the 
total power variability for every object in the MWA field 
of view at 2 s cadence. Very strong detections were trig¬ 
gered a t the position of PKS B 2322—275 (a redshift 1.27 


quasar; McCarthy et al. 1996) during a two minute pe¬ 
riod (13:30:00-13:32:00 UT). This variability prompted 
the detailed follow-up analysis described below. 

2.1. Real-Time System analysis 

Once the initial variability from PKS B2322—275 was 
identified, we searched for variability from other objects 
in the field. We an alyzed data from the MWA Real- 
Time System (RTS; |Mitchell et al.||2008| |Tingay et al. 
2013), which derives calibration solutions including flux 
densities for hundreds of sources over the MWA field-of- 
view during each integration. Eor each of the sources 
measured by the RTS (down to 155-MHz flux densities 
of about IJy/beam, depending on the position in the 
MWA primary beam) we examined the modulation dur¬ 
ing the observation. Most of the sources had no sub¬ 
stantial exc ess variability (confirming the results of |BdI| 
et al.|j2014| and |Hur ley-Walker et al.||2014|). However, we 
did identify a second source that had significant (> 30% 
over a 2-min interval) modulation, which we identified as 
PKS B2318—195. We focus on these two sources; both 
are bright and have no nearby neighbors to corrupt the 
flux density measurements (Eig. [^. 


2.2. Imaging Analysis 

To refine our analysis, we retrieved the observations 
from 2014 November 06 11:30 UT to 2014 November 06 
14:14 UT from the MWA archive. The processing fol¬ 


lowed stand ard MWA procedures (e.g.. Hurley-Walker 
et al.||2014|). We calibrated the data using an observa- 
tion of 3C 444, then rotated the phases of the visibilities 
to center the field on PKS B2322—275. We imaged in¬ 
dividual 2 s integrations in th e XX and YY instrum ental 
polarizations using WSClean (Offringa et al. |2014), per¬ 
forming 200 CLEAN iterations on each and making images 
with 3072 X 3072 0.'45 pixels, with synthesized beam full- 
width at half maximum of 2.'3. Einally, we corrected 
the instrumental polarization to Stokes I (total inten - 
sity) using the primary beam from Sutinjo et al. (2015). 
Eor a single 2-s integratiom the average image noise was 
0.16Jy/beam. In EigurefPwe show a 2 min integration 
with a slightly wider field-of-view, which has a rms noise 
of about 0.05 Jy/beam 


http://www.mwatelescope.org/astronomers 


We then used Aegean (Hancock et al. 2012) to mea¬ 
sure the flux densities and positions of PKS B2322—275, 
PKS B2318—195, and a number of other sources with 
comparable (4-10 Jy) flux densities. Both targets had 
mean flux densities consistent with catalog values (5.8 Jy 
at 160 MHz for PKS B2322 -275 and 4.8 Jy at 155 MHz 
for P KS B2318—195; |Slee||1995l [Hurley- Walker et al.| 
2014|. Eor each integration we determined the mean 
fractional variation — caused by residual instrumen¬ 
tal/ionospheric effects — of a number of other sources 
and divided it from each integration. We plot the result¬ 
ing time series in Eigure|^ We also show the modula¬ 
tion index (rms divided by mean over a 2 min interval) 
for PKS B2322-275 and PKS B2318-195. Eor com¬ 
parison, we plot the flux densities of one of the refer¬ 
ence sources (PKS B2331—312, also shown in Eig. 
and of PKS B2322—275 on another night with a com¬ 
parable data-set (2014 October 16, chosen randomly). 
During the quiescent periods both PKS B2322—275 
and PKS B2318—195 have modulation comparable to 
that for the reference source and PKS B2322—275 on 
the comparison night (roughly 5%, consistent with the 
measured noise properties). But at the peak, both 
PKS B2322-275 and PKS B2318-195 show modula¬ 
tion of > 100% peak-to-peak on 2-s timescales, cor¬ 
responding to rms variations of > 25% averaged over 
2-min. Eor PKS B2322—275 we see that the bright¬ 
est spike occurs first at the upper edge of the band¬ 
pass before progressing down, giving large apparent 
swings in the instantaneous spectral index. The varia¬ 
tions for PKS B2322—275 appear roughly 20 min before 
those in PKS B2318—195, which is about 9° north of 
PKS B2322—275 and at a slightly higher solar elonga¬ 
tion (e = 114° for PKS B2322—275 vs. e = 118° for 
PKS B2318—195) appears to start 20-30 min after that 
for PKS B2322—275. There is no indication of increases 
in the amplitude of position fluctuations during the pe¬ 
riod of most intense variability, with fluctuations of about 
5". 

In Eigure we show temporal power spectra for 
PKS B2322-275 and PKS B2318-195 (again with 
PKS B2331-312 and with PKS B2322-275 from an¬ 
other night for comparison) created from the parts of 
Eigure with strong variability. Both PKS B2322—275 
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Fig. 1.— MWA 155 MHz image of the EoRO field, centered on PKS B2322—275. The image is a single 2-min integration from 2014 Novem¬ 
ber 06 and is roughly 17° on a side. We show 1 deg^ boxes around PKS B2322—275, PKS B2318—195, and PKS B2331—312, with expanded 
views around these sources in the insets to the right. 


and PKS B2318—195 have considerably more power at 
frequencies above roughly 0.06 Hz than the comparison 
observation. Below this point the power spectra are sim¬ 
ilar, suggesting that this is the transition to “normal” 
variability (likely dominated by the ionosphere, which 
should affect all sources roughly equally). Above the 
transition both scintillating sources are similar, declin¬ 
ing with spectral slopes of roughly —1. 

3. DISCUSSION 


3.1. Origin of Variability 

We can reject intrinsic variability since it appears for 
more than one source at roughly the same time. More¬ 
over, light travel time arguments for an AGN rule out the 
2-s timescale we see, and the extreme apparent spectral 
swings do not resemble any known source behavior. In¬ 
strumental effects are also ruled out, since the timescales 
for instrumental changes (30 min) are far longer than the 
fluctuations we see and only two sources vary. We are 
then left with propagation effects, i.e., scintillation due to 
inhomogeneous, turbulent plasma along the line of sight. 
We must consider whether scintillation is of interstellar, 
interplanetary, or ionospheric origin. 

Refractive interstellar scintillation causes much slower 
(minutes to months) tens-of-percent variations in the flux 
densities of compact radio sources , with the tim escale 


increasing at lower frequencies (e.g., Rickett 


can be clearly ruled out. Diffractive scintil 


1986), hence 


ation can be 


substantially faster ( Cordes fc Rickett]|1998 ), but would 


require sizes that imply brightne ss~temperatures of ^ 
10^^ K — impossible for an AGN (Readhead 1994). 

Ionospher ic scintillation a ffects sources with angular 
sizes < 10' ( Loi et al.|2Q15b ), causing joint position- and 


amplitude ffuctn ations on timescales of > 20 s (Thomp¬ 


son et al. 2001), with the former a significant source 


of position wander in M WA data (Hurley-Walker et al. 


2014[ Loi et al. 2015b) and the latter usua lly not very 


significant (although see Fallows et al. [Mil ) Therefore 
ionospheric refraction should affect the thousands of un¬ 
resolved sources in each image. In contrast, here the 
very strong amplitude fluctuations without any position 
shifts for unresolved sources and the rather small number 
of scintillating sources argue strongly against scattering 
by ionospheric plasma. 



PKS B2318—195 are often seen to scintillate with Ooty, 
consistent with angular sizes of < O'.'l (at 327MHz), 
while Ooty sees very little (if any) scintillation toward 
PKS B2331—312. Note that Ooty does not observe past 
6 = 100°, so it could not monitor PKS B2322—275 during 
this event, although it did observe a number of nearby 
sources at e > 90° (Figs. and [3 that showed enhanced 
scintillation on 2014 November T)6 at frequencies > 2 Hz. 

Moreover, the two scinti llating sources are d etected at 


frequencies up to 20 GHz ( 

Murphy et al. 

2010|), suggest- 

ing the presence of compact structurd 

pu 

(Ghhetri et al. 


density at 155 MHz (Hurley-Walker et al. 2014) but is not 


PKS B2322—275 is one of only two sources in this field 
brighter than 5Jy at 180 MHz which is unresolved at 20 GHz 
(determined by its 6km-baseline fiux density). PKS B2318—195 
is partially resolved at 20 GHz, suggesting some arcsecond-scale 
structure but also significant compact emission. 
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Fig. 2.— Bottom section: flux densities of PKS B2322—275 on 2014 November 06 (blue), PKS B2322—275 on 2014 October 06 (black), 
PKS B2318—195 on 2014 November 06 (red), and PKS B2331—312 on 2014 November 06 (green) measured with the MWA. The flux 
densities were measured from individual 2 s integrations, with uncertainty of about 0.2 Jy. Note that the y-axis ranges for all four panels 
are the same. Bottom panel: the modulation index (rms flux density variation over a 2 min interval divided by the average flux density for 
that interval) for the same data-sets. The expected uncertainty on the modulation index under the assumption of no extra variability is 
< 1%. The top section shows the flux density of PKS B2322—275 along with the dynamic spectrum, measuring the flux density across 8 
different sub-bands (3.84 MHz each, with frequency increasing from the bottom) for the period indicated. The gaps are from the intervals 
between subsequent 2-min observations. 


detected above 5 GHz (implying a steeper spectral index 
a ~ —1.2, with Sy (X and hence likely has more ex¬ 
tended emission. This can explain why PKS B2322—275 
and PKS B2318—195 were observed to be the only sig¬ 
nificant bright scintillators. That would not be the case 
for an ionospheric origin. 

Finally, NASA’s Omniweb databas^^ suggests that 
only weak auroral and ionospheric activity occurred on 
both nights, with no clear CME arriving at Earth, mod¬ 
est Krp indice^^ (global average magnetic activity) of 
about 2, and temporal behavior over the whole field typ- 


(S. T. Loi 2015, 
le explanation for 


ical of quiet ionospheric conditionj^ ( 
pers. comm.). Absent any other viable 
what we see, and since what we see agrees with all of 
the known properties ( amplitude, time-sca le, and source 
selection) of IPS (as in|Hewish et al.||1964), we conclude 
IPS the most likely origin of the anomalous variability. 
This is largely confirmed by the Ooty Radio Telescope 
observations. 


3.2. Variability Properties 


See http://omniweb.gsfc.nasa.gov/ 

See http://www.ngdc.noaa.gov/stp/GEOMAG/kp_ap.html 
23 While we do see a modest enhancement in the auroral ac¬ 
tivity (given by the AE index) that peaks near 13:25 UT, the en¬ 
hancement is not very notable, there are multiple enhancements 
that day, and there are even more signiflcant enhancements on our 
comparison night which are not coincident with scintillation. 


The modulation of PKS B2318—195 appears to start 
20-30 min after that for PKS B2322—275, and that the 
period of enhanced modulation lasts for about 30 min 
for both sources (although some of the modulation may 
extend past the end of our observing period). Assum¬ 
ing the modulation is caused by IPS from a CME, the 
elongation difference implies a speed of ^ 500kms~^ 
for a distance of 0.1 AU. This speed is consistent with 
those inferred from Ooty observations on 2014 Novem¬ 
ber 06, where a number of sources at elongations > 90° 
that were in the same direction as PKS B2322—275 and 
PKS B2318—195 showed significantly enhanced scintilla¬ 
tion (see Eig.j^for an example). Eor a solar wind speed of 
500 7;5ookm s~ , we infer a minimum transverse size of the 
CME of 9 X 10^7;5oo km from the duration of the modu- 
latio n. However, given typical radial expansion of CMEs 
(e.g., Manoharan et al. 2000) we expect radial sizes of 
0.2-0.4AU. 'Therefore we believe that the duration of the 
IPS is set by the thickness of some sub-structure within 
the putative CME. 

As there is significant variability down to our sam¬ 
pling of 2 s, we infer CME structure on scales down to 
<10007;5 oo km. We see an apparent transition from a flat 
power spectrum to a decline in Eigurej^ This may be as¬ 
sociated with the Eresnel scale, above which IPS is sup¬ 
pressed (Manoharan & Ananthakrishnan| 1990). How¬ 
ever, actually measuring the Fresnel scale and directly 
constraining the turbulent properties generally requires 
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Fig. 3.— Left: temporal power spectra of the flux densities for PKS B2322—275 (blue, with 62 min of data), PKS B2318—195 (red, with 
37min of data), and PKS B2331—312 (green, with 235min of data), using the data from Figure[^when large-amplitude scintillations are 
evident. The power spectra are estimated using Welch’s method, averaging spectra with 256 data-points and 50% overlap. We also show 
the power spectrum of PKS B2322—275 from our comparison night (black). The cyan dashed lines show power-laws with slopes a = 0, 
— 1, and —2 as labeled. Right: temporal power spectra of the flux density for PKS B2211—388, measured by the Ooty radio telescope at 
327-MHz. The power spectrum for 2014 November 05 (blue) shows substantially less power above 1 Hz than that for 2014 November 06 
(green; note that both power spectra have been normalized below 0.1 Hz). The model fit (red dashed line) is for a solar wind velocity of 
PS 500kms“^. 


considerably faster sampling. Instead, since the rms fluc¬ 
tuations are < 100%, it is likely that the sources are par¬ 
tially resolved by IPS, the change in spectral slope could 
have more to do more with intrinsic source sub-structure 
(this would be consistent with the change in scintillation 
with radio frequency in Fig. |^. 

3.3. Identification of Candidate CMEs 

At the time of the observation, PKS B2322—275 was 
at e = 114° and a position angle of about 114° (east of 
north). So any plasma structure that caused the scin¬ 
tillation would have to have been outside the orbit of 
the Earth and below the ecliptic (Fig. |^. With typ¬ 
ical speeds of 500-1000 km s~^, we searched for CMEs 
that were observed to occur around 3-5 days before 
our observations. Among the CMEs identified by Com¬ 
puter Aided CME Tracking software (CACT us; |Rob- 
brecht, Berghmans, & Van der Linden 2009) we iden¬ 
tified a number with reasonable ejection times and an¬ 
gles. The scintillation image produced by the Ooty 
Radio Telescope (Fig. for 2014 November 05 up 
to 17:00 UT shows little scintillation in the direction 
towards PKS B2322—275, while the scintillation in¬ 
creased significantly after 20:00 UT and contiuing on 
2014 November 06. Without more detailed tracking or 
simulation of heliospheric structures, especially beyond 
Earth, we cannot conclusively determine whether a spe¬ 
cific CME or region resulting from multiple CME inter¬ 
actions caused the observed scintillations. 


4. CONCLUSIONS 

Plans for observations of IPS with the new generation 
of low-frequency radio telescopes, namely LOFAR and 
the MWA, concentrate on dedicated high-time-resolution 
(< 10 ms) observations of pre-selected sources. The tem¬ 
poral and frequency dependence of IPS flux variations 
help probe the stru cture of the solar wind (including pos¬ 
sible CMEs) (e.g., Oberoi fc B enkevitch||2QlQ| |Fallows 
et al.|2013 Bowman et aH2U13p. These observations are 
similar in strategy to dedicated prog rams with existin g 
facilities like the Ooty radio telescope (Manoharan 2010), 


-2 


which measures scintillation of 400-800 sources each day 
distributed over e < 100°. 

Assuming the variability we see is related to IPS, 
we have discussed a different regime. First, we used 
commensal night-time observations with the MWA, so 
e > 90°. Second, the sources observed at facilities like 
Ooty correspond to 155 MHz flux densities 5'i55 > 10 Jy. 
which have an areal d ensity of 100 sr“^ = 0.03 deg 
(Williams et al. 2012). Over the 600deg^ MWA field- 
of-view that means we would detect at least 20 of these 
bright sources. In fact, we can sensitively probe flux 
density variations of 20% in a single 2-s integration for 
sources as weak as 1 Jy, which have a sky density 30 x 
higher. Therefore, the density of sight-lines on the sky 
available to probe IPS is a factor of five finer, and in fact 
weaker sources are observed to vary (at lower confidence) 
along with the sources discussed here. While the 2 s inte¬ 
grations used here average over much of the higher-time- 
resolution IPS struct ure, quasi-real time analysis (as in 


Mitchell et al. 2008) can be used to identify scintillat 
ing sources which can then be followed at higher time 


resolu tion using dedicated observations (Tremblay et al 


2015). 

One difficulty in interpreting our results, though, 
comes from the fact that we are looking outside the 
Earth’s orbit, where the geometry does not allow easy 
identification of the distan ces to turbulent structures (cf., 
Armstrong & Coles 1978|). Even so, we can still exam- 


ine the statistics of IPb-like variability over multi-year 
timescales, correlating it with dedicated ground- and 
space-based solar observing. We can also probe the frac¬ 
tion of compact (<arcsec) emission from the roughly 10^ 
bright, high-latitude (where s catter broadening will not 
matter; Cordes & Lazio |2002 ) sources to be observe d in 
various MWA surveys (cf. hundreds of sources in Lit 


tie fc Hewish||19681 jeyaknmar et all2QQQ[ although see 
Hajivassiliou & Duffett-Smith 199Q[ ), providing a useful 


adjunct to calibrator surveys and allowing connections 
between spectral and spatial properties of sources. 

Finally, we note that the existence of strong IPS in 
nighttime observations of a primary field for the MWA 
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2014 Nov 05 07 UT-17 UT 


2014 Nov 05 20 UT-2014 Nov 06 18:30 UT 



Fig. 4.— Ooty radio telescope 327-MHz scintillation images. Left: image for the period 2014 November 05 07 UT to 17 UT. Right: image 
for the period 2014 November 05 20 UT to 2014 November 06 18:30 UT (whic h spans our scintillation event at 13:30UT). The images show 
the level of scintillation measured by the g parameter (e.g., |Manoharan|2010|), where = 1 is a normal level of turbulence. Radial distance 
from the center shows the derived heliocentric distance of the scattering medium: concentric circles with radii of 50, 100, 150, and 200 Rq 
are plotted and labeled. North is up, and East to the left. The arrows show the directions toward PKS B2322—275 and PKS B2318—195 
starting at a distance of 1 AU (215R©). The image from 05 November shows very little IPS, while the image from 06 November shows a 
significant increase in IPS toward the South East where the MWA data are. 


EoR experiment using nominal imaging parame ters may 
be concerning . However, a detailed analysis by [Trottfc 


Tingayl (|2015|) shows that under all but extreme condi- 

ifl 


tions, irb will not affect the detectability of the EoR 
signal. 
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